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ABSTRACT 


An analytical model for replenishment at sea is formu- 
lated for two supply ships and L combatants using queueing 
theory concepts and a random walk model in the plane. Ex- 
ponential distributions are assumed for replenishment times, 
and, given the initial number of combatants to be replenished 
by each supply ship, the distribution for total time to com- 
plete the finite operation is obtained in terms of Laplace 
transforms. All possible sequences for finishing the re- 
plenishments of the combatants have been considered in the 
model, and the techniques which were developed to count the 
number of sequence possibilities are presented as an appen- 
dix. Although this model involves only two supply ships, 
it is believed that the methods used may be generalized for 


application to more complicated models. 


TABLE OF CONTENTS 
Section 
I. Introdugeion 
lize Background 
III. Description of a Replenishment Operation 
IV. Formulation of the Problem 
V. The Model 
A. Graphical Representation 
B. Distribution Development 
C. Total Time Probability Statement 
Vi. The Solmeren for Totalmiime 
A. Path Probabilities 
B. @ransterm of the Disgri bution 
C. Examples 
BIBLIOGRAPHY 
APPENDIX 
A. “Conditional Time Distrvpution for .soters 


B. Counting Techniques for Bounded Step 
Functions 


Page 


27 


oe 





Cistl? OF TVEbUSTRATIONS 
Figure 
1. Initial arrangement and flow sequence for UNREP 


2. Graphical representation for UNREP model with 
restrictive boundaries 


3. Path illustration of reflection principle about 
lower right boundary on graph 


4. Illustration of step function that hits lower 
right boundary of graph exactly once 


Page 
Js 


JL 


28 


28 


ACKNOWLEDGEMENT 


The writer greatly appreciates the invaluable assist- 
ance and guidance of Professor Paul R. Milch during the 


course of this investigation. 


ieee Ine rOauct ion 

The capability of the United States Navy to maintain 
fleets at sea, fully ready to carry out any assigned task, is 
an important asset in the ever increasing global responsibil- 
Zties Of thigugmarion. IO acmmeve this Capability, seme 
necessary logistic support for naval combatant forces is pro- 
vided by replenishment at sea. This underway replenishment 
(UNREP), however, must be accomplished without interfering 
with the primary mission of the supported force. Consequent- 
ly, the principle aim of underway replenishment is the safe 
delivery of a maximum amount of supplies in a minimum of 
time. 

An UNREP is accomplished primarily by means of intership 
horizontal transfers via rigs connecting the supply ships and 
the supported units. Normally, there is a supported unit on 
each side of the supply vessel. This necessity for working 
at close quarters makes maneuvering a critical operation. 
Speed and course changes are restricted, and speed is 
necessarily slower than normal. The increased vulnerability 
of forces while replenishing and the increased hazards 
associated with close operations further enhance the impor- 
tance of minimizing total replenishment time. 

Extensive training, revised delivery techniques, and 
newer types of supply ships are among the ways the Navy is 
presently combatting the time problem. However, the overall 
efficiency and effectiveness of an UNREP is usually directly 
proportional to the thoroughness of prior planning. A 


detailed knowledge of the limitations, capabilities, and 
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requirements of all units involved in an UNREP is essential 
to properly schedule ships for a successful replenishment. 

Aids for efficient planning and models for studying the 
effects of this planning are limited. McCullough (5) made 
an analytic approximation of the replenishment process by 
using a multi-stage cyclic-queueing model. The model con- 
Sidered M supply ships, the stages, placed in series. These 
ships serviced N combatant units, each of which passed by 
the supply ships in succession. An infinite queue with cycles 
was then assumed by allowing the combatant ships to repeat 
the operation indefinitely, and the long run (or "steady 
state") behavior was studied. The solution to this model 
provided an upper bound for a computer simulation also con- 
sidered. 

Although McCullough's model gives some insight into the 
UNREP process, it does not really represent the actual opera- 
tion of replenishment at sea. An UNREP is not cyclic in 
nature, because the sequence of Operation is not repeated. 

Tt is a finite operation with series queues in parallel. The 
number of combatant ships is fixed, and there are no new 
arrivals. Consequently, a time dependent (not "steady state") 
solution is required. 

Gordon and Copes (2) developed a deterministic model for 
the planning of a replenishment operation by treating the 
UNREP as a special case of a job-scheduling problem. The 
service times of each ship were assumed to be known, and 
general expressions were obtained for the total time to com- 


plete the UNREP and the total waiting time of the ships 
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involved. A solution was derived which considered a maximum 
replenishment force of three supply ships and only three 
combatant ships, although the techniques developed were be- 
lieved extendable to larger operations. Unfortunately, when 
Service times are predetermined and fixed, as in this model, 
the possibility of a mishap or unexpected event occurring is 
not anticipated. Rig failures, broken lines, or accidents 
could change the replenishment time for a given ship. Thus, 
it seems reasonable to assume service times are random 


vaGrables rathemeenam deterministic. 


II. Background 

The general UNREP problem is concerned with a flow of 
customers requiring services, a trait common to all queueing 
systems. However, an UNREP differs from most queueing proc- 
esses because it involves finite series queues which start 
Simultaneously at every service facility. Each facility may 
provide a different type service, and every customer nor- 
mally requires all services that are provided. The order for 
receiving services is predetermined for every customer, but 
each customer does not follow the same sequence, so that 
there is a different series queue for each customer. All 
queues are finite, and there are no new system arrivals. The 
Operation is finished when the original number of customers 
in the system has received all required services. 

There are many complications associated with the UNREP 
problem alone, but finite inputs to systems with waiting 
lines at several service facilities are not uncommon. A sim- 
ilar problem is evident at a garage when a given number of 
cars require the same services and order of service is unim- 
portant, 1-.e., Services such as tire rotations, 011 changes, 
and gasoline fill-ups. Another example materializes ina 
commercial store at closing time, when a given number of cus- 
tomers all require services at several different counters. A 
group of refreshment stands catering to a given number of 
people also falls in this category. 

The most practical examples are probably the multistage 


production processes. These processes have been studied in 
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the category of job-shop problems, and the similarity to an 
UNREP was indicated by Gordon and Copes (2). The known job- 
shop investigations to date, however, are concerned with the 
case in which all jobs must start with the first machine. 
Other constraints usually considered involve definite se- 
quences of operations and time limitations for each job. 

Bellman (1) has described a number of simple prototype 
multistage problems and touched on some of the analytical and 
computational techniques used in early investigations. Many 
of the complications involved in an UNREP were considered, 
but never were all present in the same model. 

Sisson (6) defined the job-shop process and reviewed 
several methods for sequencing in job shops. Two basic mod- 
els for the job-shop sequencing problem were discussed, but 
they were presented only as an aid for an intuitive under- 
standing of the situation. A complete solution was not found. 

In a more recent work, Smith and Dudek (7) describe an 
algorithm that yields an optimal sequence for n-jobs requir- 
ing processing through M-machines when no passing is allowed. 
A pre-scheduled sequence is assumed, times are deterministic, 
all jobs commence with the first machine, and only one series 
is considered, therby severely restricting its use in an 
UNREP model. 

Service facilities in series have also been investigated 
using queueing theory, but, although there have been numerous 
studies in recent years, the models presented cannot be sim- 


ply modified to represent the UNREP situation. Only a 
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relatively few of the studies have even considered the 
restriction of finite queues, and these studies usually do 
not limit the queue size in front of the first server. One 
of the first investigations of finite queues under these con- 
ditions was performed by Hunt (4), who derived the maximum 
possible utilization for four particular cases of service 
facilities in series: infinite storage between stages, no 
storage between stages, finite storage space between stages, 
and the case of the unpaced belt-production line. The corre- 
sponding expected number of customers in the system under the 
assumption of exponential service times was also obtained. 

Hillier and Boling (3) extended Hunt's work in terms of 
numerical results and numerical procedures that made it pos- 
sible to analyze larger systems having exponential or Erlang 
service times. But, the input process considered was again 
such that the first queue was never empty. 

Allowing only finite inputs, starting all service facil- 
ities at the same time, simultaneously allowing different 
sequences of machine usage or service, and keeping all facil- 
ities occupied, as desired in an UNREP, create difficult 
twists to the job-shop problem or series queue situation. 
Unfortunately, there are no published works, to the author's 


knowledge, which have considered problems of this nature. 
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III. Description of a Replenishment Operation 

In an actual underway replenishment, there can be several 
different types of supply ships (servers) and different types 
of combatants (customers). A complication that is immediately 
evident is that the quantity of each type of supply ship (and/ 
or the combatant) can be different. For instance, a typical 
underway replenishment group could consist of three AOs 
(oilers), one AK (cargo ship), and two AEs (ammunition ships) 
replenishing five DDs (destroyers) and one CVA (attack car- 
rier), or one AO and one AE replenishing a CVA, a CAG (guid- 
ed missile cruiser), and four DDs. This latter composition 
is typical of the frequent UNREPs in the South China Sea. 

Each supply ship is capable of servicing two ships 
simultaneously (excluding helicopter operations). The ser- 
vice rate is normally different for each type of server and 
also varies according to the type of combatant being serviced. 
The combatant order for replenishing is predesignated so that, 
technically, lines are formed behind each server when the 
operation commences, but the service sequence is not the same 
for all customers. Each combatant replenishes at one supply 
ship and then proceeds to another line, replenishing from 
each type of server only once. The operation ends when all 
combatants have been replenished. 

The many complications of the UNREP problem make it 
amenable for computer simulation, but computer programs can 
be expensive to run, and computers are not always available 


to the planning staffs. Therefore, it seems desirable to 


is 


develop an analytical model which could realistically approx- 
imate an underway replenishment operation. The model describ- 
ed in the following pages is presented as a first approach to 


the problem. 
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IV. Formulation of the Problem 

Although a model for the general case of m servers and 
n customers is desirable, for simplicity, attention is re- 
stricted to the case of two different servers and a given num- 
ber of customers. Each customer is serviced by both servers. 
The distribution of time for completion of an UNREP is ob- 
tained. Having obtained this solution, it is believed that 
the problem may be generalized to include a larger number of 
supply ships. 

It should be noted that the solution to this problem is 
still a practical one. Many UNREPs presently conducted in 
the South China Sea involve one or two supply ships, each 
replenishing any number of ships of the same type. 

The force composition being investigated is two supply 
ships, A and B, and L combatants of the same type. Each 
supply ship replenishes only one combatant at a time. When 
the UNREP begins, the L combatants are divided into two ser- 
vice lines, M ships waiting for server A and N = L - M ships 
in B's service line. No additional ships join the queues 
once the operation has started. 

The order of service is predetermined. When a combatant 
has been serviced by A, that ship joins the queue behind B. 
Likewise, a combatant proceeds to A's line when replenishment 
is completed from B. The operation is finished when all 
combatants have been serviced by both supply ships. Figure 
1 indicates the initial queues and the flow from one line to 


the other. 


nes: 


M+Ne=b 





Figure l 


Fach server acts independently of the other server, and 
customer services are independent, so that all service times 
are independent. Transit times between supply ships are 
assumed identical for all customers and considered negligible 
in comparison with service times. The two servers have ser- 
vice times exponentially distributed, with service rates of 
A and u, respectively. It is hoped that this restriction may 
be reduced in future studies. The exponential distribution 
frequently does fit many realistic queueing situations. How- 
ever, to insure a more accurate model, the service time dis- 
tributions should be determined from known operational data. 
This is a problem for further study and is not considered 


here. 
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V. The Model 
A. Graphical Representation 
Given L combatants to be replenished by both of two 
supply ships, A and B, the distribution for the total UNREP 
time can-be determined with the aid of two variables, U, and 
re represents the number of combatants that A replenishes 
by time t, and V 


V 
, denotes the number that B replenishes by 
iLme “te 

Stalrtingmart che OFlGiiieor a Graph, aeNOoLltzontal steer 
Of ott length wink be made to the right if A finishes re- 
plenishing a combatant before B; if B finishes first, a verti- 
cal step of unit length will be made upwards. The second, and 
subsequent, steps will also be to the right or up, depending 
on whether the next customer finishing is serviced by A or B. 
With the abscissa of the graph as UL and the ordinate as Ver 


the coordinates (U V,) indicate the state of the UNREP at 


tC 


VE 
L=N+M (Aes) 
— 
5 “a 
0 M | >Me*Nn 


Figure 2 
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any time t, and the UNREP process may now be regarded as a 
random walk in this plane. 

A specific sequence of steps from the origin to (L,L) 
will represent an entire UNREP process. However, not all 


paths are possible. If A has replenished the M combatants 


initially waiting in his line, UL = M, and B has not com- 
pleted any replenishments, ss 0, then A will be idle until 
B has replenished one ship. Likewise, if Ue 1 and 
V, = 1, A will again be idle until V,_, = 2. This leads to the 


ce e 


restriction that Us = M+ V and, since the total number of 


to 


replenishments by fais Li, Ul hae < LL. Similan 


ie 


Vy SN UL S&L. Figure 2 dépicts thé graph with eee cies 


fe 


tive boundaries. 


B. Distribution Development 
Now examine a specific UNREP sequence from start to 
finish. Let X represent the service time of A and Y the 


service time of B. By definition, let 


P(X and 
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When the process begins, the first step on the graph 
in the UNREP sequence will be either horizontally to the right 
if A's service time, X, is less than B's, Y, or vertically up 
if the situation is reversed, Y is less than X. It is shown 
in appendix A that the conditional distribution of service 
time for A, given X is less than Y, is again exponential, but 
with parameter A + uu, 1.e., 


=. G7 Mees 


P(X = Xs aye eee eee A 
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Likewise, the conditional service time distribution for B, 
given Y is less than X, is also exponential with parameter 
XA + uw. Therefore, for the first step, 


p(y <tly <x) = P(x <t|[x<¥y) =1-e ATH Y 


|v 
© 
® 


The second step of the sequence will be examined 
from the point (U, —r, Vy = 0). Again, the choice for this 
step is either to the right or up, depending on whether A or 
B finishes first. Since A iS Servicing a new customer, the 
distribution of service time for A is known. Bis still re- 
plenishing his first customer. However, by employing the 
memoriless property, P(X < T|X < t) = P(X < T-t) when t < T, 
of the exponential, the distribution of the service time that 
is left for B remains exponential with rate u. Therefore, 
the problem is exactly as before, and the conditional distri- 
bution of time for the second step 1S again exponential with 
parameter A + U, given that it is a vertical, or given that 
it is a horizontal, step. 

In general, this conditional distribution is true 
at any point except those points on the boundary, i1.e., when 

. = SS = eo «tf eee Ve 
then the path position for the process is on the right bound- 
ary, and B is working while A is idle. Therefore, the next 
step must be up, and the distribution of time is the same as 
the distribution for Y, exponential with parameter tu. This 
is again due to the memoriless property of exponential dis- 
tributions. Similarly, when Mee = Ne U, eee ene location 


is on the upper boundary, and the next step must be to the 


right with the same distribution of time as X. 
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Let al denote the event that a given route has i 
positions on the upper boundary and j positions on the right 
boundary. Then, when ae occurs, 2L - (i + j) positions 
along the path are not on a boundary. 

The total time T for a given UNREP sequence is, 
therefore, just the sum of the times between each of the 
positions on the path. Let Z be the random variable with 


Bee OA 7 ut, 


P(Z ae) eee then 


I = - 2 FR. + Y,, DEsevent Ee. ecemma. 
FTG 5) ah 5 1) 


The subscripts on Z, X, and Y indicate the number of times 
each variable is summed, l.e., X. is the sum of i exponential 
random variables, each with parameter i. 
C. Total Time Probability Statement 

Unfortunately, a given Pa Can occur in many ways, 
and, in addition, there are many combinations of 1 and j that 
are also possible. However, once the probability of event 
Be is determined, the unconditional distribution of UNREP 
time can be ascertained by applying the theorem of total 
probability. The unconditional probability statement is, 


therefore: 


Paien £) = AP. < tlE.) * P(E;.) 
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View (Che Soluticonmto, Total Time 
A. Path Probabilities 
Consider once more the step sequence representing an 
UNREP process. Starting at the origin, the first step in the 
sequence will be to the right with probability p, or up with 
probability 1 - p= q. It can be shown (See appendix A.) that 
r 
Fi 


LS 2 ar 


23s 
u + 


and gq = 7 Tie Then, considering any point along 


the route except those points on the boundaries, and employ- 
ing the memoriless property of the exponential distribution, 
the probability that a step will be to the right is p and 
that a step will be up is q. 

When a position is on the right boundary, the next 
step must be up with probability one. Similarly, from posi- 
tions on the upper boundary, the next step must be to the 
right with probability one. Consequently, when Bi SECs ; 

i steps to the right and j steps up are made at the bound- 
aries. Since L steps are made in each direction, L - 1 
steps are to the right, each with probability p, andL - j 
vertical steps are made, each with probability q- 

For a given Boas however, several different step 
sequences are obviously possible. Therefore, for specific 1 
and j, the number of possible routes must be counted before 
the probability of E can be determined. Utilizing combina- 
torial analysis and the reflection principle, a counting 
technique is developed in Appendix B which yields Boas defined 
as the number of possible ways in which ee can occur. The 


results are aS follows (Note, ~C(n,x) = n!/r'(n-1) tie 


Zall 


00 
Page Rio 
3 os 
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Since the random variables pn 2 aire 2 


= Q 
= C(2L-i-1, L-l1) = C(2L=1=1, Mea) =se27-1-1) N-1), 
Juo = los Zaee 


= Ci 2G jig L-1L) = C(2L-j-1, Mea) = C (2a, N=j), 


= C(2L-k-1,M-k+1) + C(2L-k-1,N-k+1l) - C(2L-k-1, M-k) - 
C(2L-k-1, N-k), where k = i+ j = 2,3,..., 
max (MjN) aa ander. > sie 


Su ai the probability of Bee is determined below: 


L-i L-j 
13 ijP 4 
Transform of the Distribution 
The necessary equations and expressions for determin- 


total time distribution are now known and age sims, 


briefly below: 


pe ct Xs cts Yau given Bags 
. Ce 
pa) Je . i 
—o) —— PE us. ee ae Ol hee ma os, 
a) ib ( ee ( Tags | i? 


Muon are mutually 


independent, the Laplace transform, denoted f*(s), will be 


used to 


fr (s) = 


determine the total time distribution. 


) = () eta (s) 
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iis) 


Oem 0 5-0) ess) (ura) 
~rs _ wU+ts | 
Let a = eS and 8 = oR Then substituting 


the values of oe and evaluating, if A # u: 


LL i me - 
Enis) = = ° ) C(2L-j-i, L-1) {a J 48 3] 
(A+u+s) j=1 = 


. -j -u+s =a hype * 
= ees 7-1, N= aoe, eae: —-- =. | 





j=1 \E Ses) (u-d) 
- C(2L-j-2, N=3) [a 97? AES) = sot ts) | 
M | 2 
: “Sy ea =~ ieee SD) 
fev j-1, M i) fa Gap + © O-aypqeay 
- C(2L-3-2, m3) a 9 # AES) z tts) | 
IETS) tc (2L-2, Nom) Gi on—2.) Moin 


L 
fF*¥(s) = aorta) ™| ) C(2L-j-1, Ee, ~ 
jee 


ee (2L=j-teeN—4) 4G ( 2-7 ae “ay 
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e r IG ° o4 0 = -1i 
“ layrs) [hc -2 M-j) jo? 
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+ ) C(2L-3-2, N-3)ja 7 
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T5 


C(2L-j-1, N-j)ja ? 
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For a given number of combatants, L, the effects of 
changing M and N on total replenishment time can be determined 
by differentiating fr (s) and setting s = 0. This yields the 
negative of the expected replenishment times for easy compar- 
ison. 

C. Examples 
1. Assume the simple case of two servers and two 


customers with M ="Ne=—1. Then) see, 


eee 2 a: 7, 
fA(s) = —~~—74 J €(3-4,1)fo I+8 7) - a -[C(2,0)+C(2,0)] 
(Atuts) © (j=l 
i Z 
a aba = +s) 
2 f feces srfa etsy + Reape 
Be (Qos nee e ee (OE eer (se 
C(2-j,1 300 Ges: uss 
= Nae : (u+A+2s)- = | yA (utAt+2s) ]’ 
(Atuts) 7 (ats) (u+s) 7 (Atuts) (At+s) (uts) 
Differentiating and seubing Ss = 507 
E(t) = 20) = 


2. bet L = 10 with M= N = 52 Then weer eee 


0 


One -4. 973 = 
£*(s) = —SE 6) J c(19-5,9)(a 348-7) - 20° "c(18,4) 
(A+ut+s) j=1 
- 2 ; €(19-4,5-4)-|a_- Ets )+ Jel ile (ats) “ )| 
j=1 SL Ges Leer Oem. 


+ o(18-5,5-5) [a 72 (ASS) - Thats] 
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If A = 2 per hour and u = 4 per hour, after differentiating 
and setting s = QO, 


(Tl) =eet Siete 5.15 Hoursee 


2D 
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APPENDIX 
A. Conditional Mime Distrib@eron for Steps 
The random variables X and Y have exponential distribu- 
tions with parameters X and u, respectively. The conditional 
distribution of X, given X < Y, is developed below: 


P(t oo) 


om dk ly =i “Sa < 


l 
rg 
. 
V 
ct 

4 

A 
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P(t < X < Y) f(x) fy(y) dx dy 


1 
N—- 





clei “ 7 , t > 0 
liveeneo, then 
P(X < Y) = x 7 
Therefore, P(X et) x <  V)— ea 35 ude ed) 


B. Counting Techniques for Bounded Step Functions 

Step functions having unit steps which are vertically up 
or horizontally to the right are the only functions permitted 
in a graph with the following boundaries: x=M+t+ty, x=L, 
y = N+ x, and y = L. The boundary x = M+ y will be the 
lower right boundary, and the combined boundaries 
xX =M + y"= Lo will be referred EoO.as Ehbewerght boundary ime 
upper left boundary is y = N + x, and the upper boundary is 


the combined boundaries y= x +N < L. 


Z7 






a (y+M, x-M) f 
Original /-* 


Path 
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Ga 2 M L=M+N 
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Path illustration of reflection principle about x=y+M,. 







» Lee) 

Path after mandatory 
step deleted and 
subsequent steps 
deeuped one unit. 


M MH 
Figure l 
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AP) 


Illustration of step function that 
hits lower right boundary once. 
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First a technique for counting the number of step 
EUNCETONS Lrem 00h) to (x,y) < x < Land N Sey Se, 
that do not touch the upper left or lower right boundaries 
must be developed. To do this, the number of paths or routes 
(both terms used interchangeably with step functions) that 
have a point in common with these boundaries will be sub- 
tracted from the number of routes that are possible without 
considering any boundaries. 

The number of steps required to go from (0,0) to (x,y) 
by any route is x steps to the right and y steps up. The 
number of possible combinations of x + y steps with x steps 
femthe, riqnagesceGtze+, y ~ooe= {xt /xivin This is also 
the number of routes that are possible without considering 
any boundaries. 

Next, consider a path that hits the lower right boundary 
of the graph as exemplified in figure 3. If the part of this 
reute from the first point of Contact with the boundary te 
(x,y) 1s reflected symmetrically about the boundary, then the 
reflected route (dotted line in diagram) has the same number 
of steps to the point (y + M, x - M) as does the original 
step function to (x,y). AlS@7eaf any Ppathewenalr, ome ececcesuce 
(y + M, x + N) is reflected symmetricaliy from the first 
point of contact with x = y + M, the reflected route proceeds 
to (x,y). Therefore, the two sets of routes are equivalent, 
and it is only necessary to count the step functions from 
(0,0) to (y + M, x - M) in order to determine the number of 
routes touching the lower right boundary. This number is 


C(xt+y,y+M). 


2g 


Similarly, by using the reflection principle about 
y = x + n, the number of routes touching the upper left 
boundary is C(x + y, x + N). Therefore, the total possible 
number routes from (0,0) qgto (xgy), Meg x < Upset 


that do not touch either of the boundaries is 
C(x + y, X) - C(x + y, y + M)>= Gixs+ y, x Fe 


Note that this relation is dependent on M and N, the inter- 
section of the boundaries with the axes, so that the bound- 
arles may be shifted and the total number of routes not 
touching the new boundaries may be determined by the above 
method. 

It can be shown that, in general, the number of possible 
routes from any point (XO, ve) to(x,y) that do not touch the 


boundaries xX = y + Mand y = x + N is 
C(x-xX oty-y. 1 X-X)) = C(X=X TY Vo eee - C(x-x ty-yo  X-YStN) . 


Now define Ris as the number of possible paths from 
(0,0) to (L,L) that touch positions on the upper boundary i 
times and the right boundary j times. Since a diagonal step 
is not permitted, every path must pass through (L - l, L) or 


(L, L -.1) to get to (L, Lb), and so ke is obviously zero. 


00 


Consider Boys next, where j = l, 2, ..., L. In thisgeascc ae 


routes hit the upper boundary. 


To formulate a method for finding R in general, it 


05 
will be shown that a set A, consisting of routes hitting 
positions on the right boundary exactly j times, is equiv- 


alent to a set B, consisting of routes that proceed to 
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(lL —"f))L - J) witheoutehitting the boundaries = y + m + 

j Sete seL andey =e N < Lb. This being true, the technique 
developed to count routes within specific boundaries can be 
applied directly. 

If j = 1, set A paths pass through (L - 1, L - 1) with- 
out hitting either of the original boundaries and then move 
tothe srtghtaeecee Bo paths alsemproceed to (Lb -— 1, Ley 
without ha tting@ene boundaries ea yet M =< E and 
y=x+N< L. The two sets are obviously equivalent be- 


cause they are the same sets, and 


mil hy al 2 re aS al OA rt LN) 


If 3 is greater than one, however, the equivalence is 
not so obvious because j different types of paths are pos- 
Sible. For descriptive ease, "type K" will indicate the 
type of path that kits K positions on the boundary x = L and 
j - K positions on the lower right boundary. For example, 
let j = 2, and consider set A. Two types of routes are pos- 
Sible, type I and type II. A type II route passes through 
(L - 1, L = 2) without hitting any boundary and proceeds to 
(L, L - 2). A type I route hits the lower right boundary 
exactly once, proceeds to (L - 1, L - 1), and then to the 
real glare 

For j = 2, since all routes in set B proceed to (L-l, 
L-2), type II routes in set A correspond to the same routes 
in set B. A type I route, however, hits the lower right 
boundary once, and the next step after hitting the boundary 
must be up. If this mandatory step up is deleted from con- 


sideration, and all subsequent steps of the path are dropped 
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one unit down, as indicated by the dotted path in figure 4, 
the path proceeds to (L - 1, L - 2) and then right, never 
touching x = y +M+t+4I1. Therefore, all routes in set A 
have corresponding routes in B. 

In set B, when j = 2, all routes proceed to (L - l, 
L - 2) and never touch x = y + M+ 1. Consider a route that 
hits x = y +M. If the route is raised one step after hit- 
ting this line for the first time by inserting an extra 
vertical step, then the route proceeds to (L - l, L - 1) 
and corresponds exactly to a type I route in set A. The 
remaining routes in B do not hit either x = y+M+l1 or 
x = y + M and still proceed to (L - 1, L - 2). These routes 
correspond to identical type II routes in A. Consequently, 


the two sets are equivalent, and, using set B, 


Roo = C(L-1+L-2, L-1) - C(2L-3, L-24+M+1) - C(2L-3, Lien) 
= C(2L-3, L-l) - C(2L-3, N-2) - C(2L-3, M-2). 
In general, when i = 0, j types of routes are possible 
among the R,. routes in set A: types I, II, ..., andj. A 


07 

type j - 1 route hits the lower right boundary one time and 
then proceeds to (L-l, L-j-1l). If the mandatory step up 
after hitting the boundary is deleted and the remaining part 
of the route is dropped one unit, the path will proceed to 
(L-1, L-j) without hitting x = y + M + 12 For age eee. 
route, I < K < j, j - K steps must be up when Ehem@lower @amime 
boundary is hit. If these mandatory steps up are deleted and 


the remainder of the route is dropped one unit following each 


deleted step, the path proceeds to (L-l, L-j) and corresponds 


az 


teomaeroute In Bo othatvaeecs not hit x = y +M%+ K. Therefore. 
all routes in A correspond to routes in B. 

In set B, consider all routes that do not hit x = y + M. 
These routes correspond exactly to the type j routes in set 
A. Next, consider a route in B that hits x = y + M + j-1-K 


and does notyhriw< — y + M + 75— Kia j - l. If a ver- 


[| A 


tical step is added to each route immediately following the 
first and subsequent times the boundary x = y + M is hit, and 
the remainder of the route is raised one unit for each step 
added, the raised route will correspond to type K route in 
set A. Therefore, all routes in B correspond to routes in A, 


and set A is equivalent to set B. Consequently, for any j, 


Ro = C(L-l+L-j, L-1) - C(2L-1-4, L-j+M+j-1) - C(2L-1-j, L-1+N) 
= C(2L-1l-j, L-1) - C(2L-1-j, N-3j) - C(2L-1-j, M-3j). 
By interchanging i and j}, M and N, it is obvious that 
yin C(2n=T=1, G-1) - C(2L-l-1, M=-i) = C(2L-1l=1, N=-i)- 


Consider now Roe when i # 0 and j # 0. Define Ba as 
the number of paths that hit the lower right boundary first 
and then proceed to (L-i, L), and Bog as the number of paths 
that hit the upper left boundary first and proceed to 
(Lb, L-))-Siienpa: 1) ~ 0 andi a = 0; Ba = Bas + Ba 

Let 1= 1. Then, £or Ria? the number of routes that hit 
the lower right boundary j times and proceed to (L-1, L-1l) 
must be found. Consider type I routes among the set A 
routes that hit thereignt boundary )°+)1 times and do net tie 


the upper boundary. These routes hit the lower right bound- 


ary j times and proceed to (L-1, L-1), which is exactly what 


33 


is being sought. Therefore, for R:! 


oy it is necessary only 


to find the number of type I routes for set A. Considering 


the previous set B, this is simply R less all those 


0, j+l 
routes that do not hit the boundaries x = y+M+ (jy + 1) - 
Z2=y + M - j = land ye "x + Newer proceeding eo 

(L-1, L-j-1). Therefore, 


Ri 


[(C(2L-1-(541) ,B=1) = C(2L=4=2)N=4-1) = C(20=jS2eee ee 


[C(L-1+L-j-1,L-1) - C(2L-j-2,L-4-1+M+ 5-1) -C (20- jae 


II 


C(2L-j-2,N-j) - C(2L-j-2,N-j-1). 

Let i be greater than one, and consider a new right boundary, 
xX=y+MseL- itil. If the number of type I routes Higa 
this boundary j+l times and not hitting y = x + N is found, 
then Ba can be determined. By the same techniques used to 
formulate on when x = y + M < L, if the right boundary is 
xeScy st Messieere > ay, “then 


RS gay = C(bratb-(j+1), L-i) - C(2L-i-j-1, L-j-1+M+ (j+1)-1) 


C(2L-i-j-1, L-i+N) 


I 


C(2L-i-j-1l, Lei) - C(2L-i-j-1, N-i-j) 
- C(2L-i-j-l1. M-j-1). 
The number of type I routes for this case is, therefore, equal 


to RO a3) less the number of routes proceeding to (L-i, 
a 


L-(j+1)) without touching the boundaries x = y + M+ (j+l) - l 
and y = x4 N. Thus) for a7 0 ae 


Re 


jy = [e(2b-i-j-1, L-i) - C(2L-i-j-1, N-i-j) - C(2L-i-j-1, 


M-j-1)] - [C(2L-i-j-1, L-i) - C(2L-i-j-1, L-j-1+M+j-1) 


- C(2L-i-j-l, L-i-N) ] 


C.(2L-i-j-1, N-i-j3+1) = Cl 2b=i-—9-19 eae 
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tf 17+) = k, thengeeere 720 and 4.0, 

ee = C(2L-k-1, N-k+l1) - C(2L-k-1, N-k). 
Similarly, by interchanging M and N and i and jj, 

Ris =IG2t=k-1, Make) = C(2L=k=-1, Mo-icr 
tmegencral,thememonre, With a9] =k, 2 7 0, | -aGe 

Hee = C(2L-k-1, N-k-1) + C(2L-k-1, M-k-1) 

- C(2L-k-1, N-k) - C(2L-k-1, M-k). 

It should be noted that, although the maximum ranges of 1 and 
j are from 0 to L, the equations hold for all values of i and 
He since, if i OF j 1S greater t#an L, the cemernaterral 


forms will be zero. 
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